Diffusion as mixing mechanism in granular materials 
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We present several numerical results on granular mixtures. In particular, we examine the ef- 
ficiency of diffusion as a mixing mechanism in these systems. The collisions are inelastic and to 
compensate the energy loss, we thermalize the grains by adding a random force. Starting with a 
segregated system, we show that uniform agitation (heating) leads to a uniform mixture of grains of 
different sizes. We define a characteristic mixing time, Tmix , and study theoretically and numerically 
its dependence on other parameters like the density. We examine a model for bidisperse systems 
for which we can calculate some physical quantities. We also examine the effect of a temperature 
gradient and demonstrate the appearance of an expected segregation. 



I. INTRODUCTION 



Granular media are notoriously difficult to mix. For a variety of reasons and under rather general conditions, they 
tend to form segregated steady statesEl. For example, segregation can occur during granular flow, the bigger particles 
moving farther than the smaller ones. Difference of size (polydispersity of the grains) or difference of material (different 
kinds of grains) produce different geometrical or physical properties. Segregation can also be due to percolation! where 
the small grains fall through the holes between the big grains leaving only the bigger particles behind. SheaitTO and 
vibration can also produce segregation. One of the best known examples of vibration segregation is perhaps the 
"Brazil nut effect" . In this case, the geometrical propertiesBo are responsable for the upward movement of big 
particles although convection processes near the boundaries can also be very importanlH_All these processes (flow, 
shearing and convection) are very common in industrial applications such as in mixersli30. For these reasons such 
mixers are efficient only for rather homogeneous materials. In the polydisperse case it is very hard to avoid segregation. 

For gases and liquids, the thermal agitation of molecules4s a natural and efficient mechanism leading to throughly 
mixed systems with homogenous equilibrium steady states!! 2 ]. We propose here, a study of a system of agitated grains 
in analogy with liquid or gas molecules at the miroscopic scale. 

Two major differences between granular materials and fluids are: (f ) the particle size compared to the mean free 
path and (2) the inelastic properties and friction responsable for energy dissipation during collisions. The question 
then is if these differences will alter the system's natural tendency to mix by diffusion. In other words, is it possible 
to use diffusion to mix grains. In spite of the dissipative collisions, it is possible to keep a granular system agitated, 
for example on an air-table or a vibrating bed. To simulate numerically such constantly agitated granular systems, 
we add an external random force to the equations of motion (see section [□]). We then analyze the grain diffusion and 
its dependence on the various parameters of the system such as grain size. We find that inspite of the dissipative 
nature of the collisions, diffusion is still a good mixing mechanism just like in fluids. 

In section || we detail the algorithm and summarize the principal dynamic equations and parameters of our system. 



We verify our procedure with the study of a monodisperse system in section [II and etablish relations which characterize 
the temporal evolution of an initially segregated system. The bidisperse case is studied in section [V. We show in 
particular the evolution of a system with ho mog eneous agitation, and also the effects of a gradient in this agitation. 
Our conclusions and discussion are in section |V[ 



II. ALGORITHM AND REVIEW 



We use an Event Driven Molecular Dynamics algorithm. The simulated .grains have the same characteristics (e.g. 
normal and tangential restitution coefficients) as measured experimentallyO. To thermalize the system, we add at 
regular time step intervals, dt, external random forces which act on every particle. There are several choices one can 
make for this force. Our choice is the following: 

Ft = m (J$/dt\ Q (2.1) 

with i = x,y (corresponding to the two directions), Q is a gaussian noise caracterized by < QQ >= Sij and m is the 
particle mass. 7/q is the control parameter which we use to increase or decrease the agitation. Experimentally, on the 
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air tableEJ, plastic disks of radius R move due to the fluctuations of the air flux acting on their surfaces. Therefore, 
since their mass is proportional to R 2 , we expect the accelerat ion to be independent of R. That is why we have 
chosen an external force proportional to the particle mass, Eq. J2.l| ). The equation of motion of a particle between 
two collisions becomes: 




In this paper Vi denotes the instantaneous velocity in the i direction and v 2 the mean square velocity. The system is 
two dimensional and is enclosed in a square box whose walls are made of grains of radius r w and are infinitely massive. 
The particle-wall collisions are taken to be elastic. Note that in these two-dimensional simulations, the particles are 
represented as spheres interacting at their equators. 



A. Macroscopic characteristics of the steady state 

The above model leads to a steady state characterized by a constant mean square velocity v 2 (t — > oo) for all 
particles. The energy loss during collisions is compensated for by the random force. For a monodisperse gas, the 
energy balance is easily calculated. The energy loss, F, per unit time in the steady state is given bytB 



r oc Lumv 



(2.3) 



where u is the frequency of col lisio ns. On the other hand, the average gain in energy due to the random force during 
dt is easily obtained from Eq. (2.2): 
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to [v 2 (t + dt) - v 2 (t)] = mr) 2 dt. 



(2.4) 



In the steady state we can thus write 



1 dv 2 

— TO — — 

2 dt 



mr) 2 . 



(2.5) 



r~2 

We can write ui ~ ^j— where I is a caracteristic length depending only on the packing fraction and on the radius of 



the grains. With this assumption and the fact that 



it — >c 



0, Eq. (2.5) gives 



v 2 (oo) oc (Ir, 2 ) 2 / 3 . 



(2.6) 



This power law is independent of the coefficients of restitution and friction if dissipation is not too large. The results 
of the kinetic theory of inelastic gases can be applied, in particular the velocity distribution can be approximated by 
a Maxwellian. The parameter rj 2 allows us to change the granular temperature, T, since v 2 oc T. Therefore, in the 
steady state, T is independent of the initial conditions. For polydisperse gases the problem becomes more complicated 
as will be seen below. 



B. Coefficient of diffusion 



Since the mean square velocity is constant, and consequently the collision frequencies too, particles have a simple 
diffusive behavior. The mean square displacement, < (r(t + to) — r(io)) 2 >, gives the coefficient of diffusion 

<{r(t + t )-r(t a )) 2 >=4Dt, (2.7) 

where to is large enough to ensure thermalization of the system. Clearly if the system is examined at a time scale 
t < 1/oj we will not observe the true diffusive behaviour. At short time, Vj is not constant due to the action of Ff and 
so the mean square displacement is not yet linear with t. For (» l/wwe verify the linear dependence of the mean 
square displacement on time, for all particles. The value of the coefficient of diffusion D found from the simulation 
is thus larger than the theoretical value predicted by a Langevin description due to the dynamics of the particle at 
short timeO. 
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C. Simulation procedure 



We study systems made of two species of grains, s and b. The radius of the particles are respectively R s and i?t>. 
The system is a square box of length L and we use the boundary conditions discussed above. The number of particles 
of each species is calculated based on the desired packing fraction C and the relative proportion of s particles x s . 



where n s and n& are repectively the number of particles s and b. 

For all mixtures, we performed two types of simulations. In the first one, the two species s and b are already 
mixed and the initial position of each particle is chosen randomly in the box by using a classical algorithm of Random 
Sequential Adsorption (R. S. A.)lia. We are careful that this algorithm does not introduce segregation in the initial 
configuration. For the second type of simulation, the two species are initially separated with the s particles on the 
right and the b particles on the left. The system is prepared such that the packing fraction is homogeneous in the 
whole system. 

In the following section we present our results in the simple case where the two species have the same mechanical 
and geometrical properties, i.e. s and b grains are of the same type. 



To test the validity of our algorithm we start with identical grains, i.e. R s = Rb- In this case, we do not expect 
segregation because the grains are identical, but we would like to verify that the thermal process is efficient and that 
particles do not agregate. In other words, after some time each grain will have visited all regions of the box. In figure 
|l| is presented the temporal evolution of the system for the two different initial configurations specified above. In 
figure [l]a the system is already mixed and in [j]b the particles are initially separated. As one can see in this figure, 
the system does not collapse and the grains are homogeneously distributed in the box. To analyse the dynamics of 
the mixture we measure the quantity N St b(t) defined as the number of collisions between s and b grains per unit time. 
The evolution of N 8 f,(t) with time gives two important results. For the initial configuration corresponding to figure 
|l|a, the quantity N St b(t) fluctuates around a mean value, 7V Si f,(oo), as seen in figure | This means that the system has 
reached a steady state in which the mean square velocity of the particles is constant and equal to v 2 (oo). Note that 
the system evolves very quickly into this steady state. We have checked that all the configurations at different times, 
t, are statistically identical and that the system remains homogeneous. There is no evidence of collapse or cluster 
formation. The second observation is that for the initial configuration corresponding to figure U>, the quantity N s ^(t) 
increases and then stabilizes at large time at the value N a ^{po) defined above. However, the mean square velocity 
of the grains, v 2 (t), reaches the steady state value v 2 (oo) much more quickly since the grains are identical. The 
knowledge of the mean square velocity is not sufficient to define the state of the system since it gives no information 
about the spatial repartition of the two species. N St b(t) is therefore the only pertinent quantity to characterize the 
homogeneity of the system. 

A few more comments about N St b(t) are in order. In the system studied above, the packing fraction and the 
velocity distribution are spatially homogeneous and constant in time (except at very short time). As a consequence 
the quantity N S: b(t) depends only on the spatial repartition of the two types of grains. The evolution of N St b(t) 
allows us to define a mixing time, T m i X . We have already mentioned that the velocity is the same for all particles and 
independent of position. This is also true for the local density. We conclude that also the frequency of collisions is 
the same for all particles and is space independent. In this monodisperse case, the dynamics are purely diffusive and 
can be characterized by a coefficient of diffusion D which is independent of horizontal spatial position x. Let us call 
8N St b(x,t) the number of collisions between the two species occuring at a position between x and x + dx at time t. 
Clearly, SN St b(x,t) is directly proportionnal to d s (x,t) and db(x,t), the densities of s and b grains at position x. The 
densities d s and db do not depend on the vertical position since the system is invariant along this direction. We will 
define do as the total density and can, thus, write do = d s (x,t) + db{x,t). do is of course independent of x because 
the system remains homogeneous. We then obtain an expression for N s> b(t): 




(2.8) 



III. MONODISPERSE CASE 



L 




(3.1) 
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The density of big particles at (x,t), d b (x,t), is described by Fick's equation, 



with the following boundary conditions: 



dd b (x,t) d 2 d b (x,t) 



dt 



d b (x) = d 
d b (x) = 
d b (x) = d /2 



We assume the solution of Eq. (3.2) has the form 



D- 



dx 2 



for < x < L/2 and t = 0, 
for L/2 < x < L and t = 0, 
for all x and t ~ oo. 



(3.2) 



(3.3) 



do 



d b (x, t) = (B m sin(\ m x) + A m cos(\ m x))exp{-\ 2 m Dt) + — 



(3.4) 



where the X m are constants. Using the conditions Eq. (3J3) gives 



d b (x,t) = ^2a k cos(\ k x)exp{~\ 2 k Dt) + 



k=0 



with 



2d (-l)- 



afe 

x _ (2fc'+l)7T 



7r(2fc+l) ' 
"'+1' 
L 



Eq. (3.1) then gives the final expression for N Sib (t), 



(3.5) 



(3.6) 



N, ib (t) cx ^ ( 1 



(2fe+ 



w exp(-2\ 2 Dt) 



^2 _ (2fc+l)^^ 



(3.7) 
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As a first approximation, we may keep only the first mode and write: 

N s>b (t) ~ iV s . b (oo) (1 - exp(-t/T„ 

I mix — 27r 2 L>' 



0), 



(3.8) 



where r m i X can be taken as the typical time for mixing. 

To check the validity of the theoretical expressions for N s>b (t) and r mlI established above, we have performed 
simulations for different values of x s and C. For a given set of parameters, we have performed five simulations 
corresponding to different initial positions and velocities of the particles for the case where the s and b grains are 
initially separated. Figure [| shows N Stb (t) versus t averaged over the 5 simulations. We see that the agreement 
between theory and numerical simulation is very good. Figure |4| shows the dependence of the mixing time, T m „, on 
the coefficient of diffusion D. To get this, we performed several simulations changing the packing fractio n a nd the 
radius of the particles in order to vary the coefficient of diffusion. Note that D was estimated using Eq. (2/7). The 
slope of the curve is exactly that predicted by the theory. As the coefficient D can be calculated from the parameters 
of the system (R s = R b , rj^, C), we can esti mate analytically and with high accurancy the mixing time. We have also 
verified the dependence of T m i X on L 2 , Eq. ( |3.8| ), and have found very good agreement too. 



IV. BIDISPERSE CASE 



We now discuss the case of a binary mixture. We will see that the size difference between the grains changes 
drastically the dynamics of the system. The grains s and b are taken of equal density and identical coefficients of 
restitution and friction, and we take R s < R b . 

We present first the case where the system is thermalyzed uniformly, i.e. F\ does not depend of the position of the 
grain. Then we will examine the case where a gradient is imposed on the agitation force. 
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A. Case of homogeneous agitation 



We will show that, in the bidisperse case, the system also evolves into a homogeneous steady state. We will see here 
that the form of the thermalization force and the initial conditions determine the evolution of the system towards the 
steady state. In the s imu lations the packing fraction is fixed to 40% and x s , which represents the relative proportion 
of small grains (Eq. fl2.S|)), is the only parameter to be varied. 



1. Evolution at short time 

Figure ^| shows the evolution of the system with time t. In the initial configuration, the two species are separated 
and the two populations have the same initial velocity distribution and therefore v 2 = v 2 . The initial local packing 
fraction, as one can see in figure [5], is the same in the whole box. 

Recall that in our simulations the surface occupied by a particle in the plane is proprotional to R 2 and its mass 
to R 3 since the particles are spherical. Since the pressure is P cx mdv 2 , where d is the density of grains, the initial 
pressure is larger for the bigger particles. The system therefore has an initial pressure difference which will govern 
its behaviour immediately after the partition is removed whereby the bigger particles, 6, compress the small ones s. 
As t increases, (see figure the density of the b particles decreases and so does its pressure. On the other hand, 
the density and pressure of the s particles increase. During this compression period we can consider the system as 
two interacting monodisperse systems. In the left part of the box ( occup ied by the larger particles), as the density 
decreases, the mean free path, I increases. We have seen in section II A| (see eq 2.6) that v 2 increases with I. The 



velocity v£ is then increasing with time. For the same reason, in the right part of the box, v s is decreasing. As a 
consequence, the pressure, which is proportional to the product of the square velocity and the density, is maintained 
almost constant in each subsystem. The pressure difference between the two subsystems remains therefore important 
and favors the compression of small particles. The packing fraction of the s grains increases up to a value around 
68%. 

It is worth noting that if the walls were inelastic, collapse would occur whereby the small grains would be squeezed 
near the wall and would loose all their energy due to dissipation. In our simulations we use elastic wall and thus 
observe a reflection of the compression wave. To illustrate this we show in figure^ the packing fraction of the s grains 
as a function of t and x. We observe a compression wave which traverses the system. On average, the small particles 
remain compressed and the big ones dilute. During this process the diffusion between the big and the small particles 
is very efficient due to the high concentration gradient. The evolution of the quantities v 2 (t) and v 2 (t) as a function 
of time is illustrated in figure 0. One can see that the velocity of small particles decreases at short time and then 
increases when the mixing process starts. At long time the mean square velocities of both species reach a constant 
value corresponding to a steady state. 



2. Mixing time 



After the compression phase, the system starts to mix. As we have done for the monodiperse case we examine the 
quantity N s j,(t). To have a good estimate of N s f,(t), we take (as in the previous section) the mean value obtained 
over 5 simulations. Figure ^ shows that we can appoximate 2V S & by N{oo)exp{— t/r m i X ). Note that the compression 
phase occurs during a short time compared to r m j x . We obtain in this way T m i x for different values of x s for R s = 0.4 
and Rb — 0.6. It appears that this time r m i x can be considered as independent of x s (see Figure ^). This mixing time 
obtained in a bidisperse system is smaller than that obtained for the same packing fraction in the monodisperse case. 



3. Steady state 

Even though the collisions are dissipative, the system does reach an out of equilibrium stationary state due to 
the random agitation force. This stationary state should be characterized by macroscopic functions which should be 
independent of time. 

Using thermalized configurations (long evolution times), we performed a geometrical analysis using Voronoi tessel- 
lation to check that no segregation exists. We calculated the number of s neighbours for a b particle and found that 
the distributions of neighbours are roughly identical to the distributions obtained from static configurations generated 
by an R.S.A. algorithm. We should point out that the distances between particles can be different from the static 
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case but the neighbourhood of a grain is the same in the dynamic and static situations. This demonstrates that there 
is no segregation. 

We now consider the distribution of the kinetic energy as a function of the radius of grains. In the case of elastic 
collisions, we can define a kinetic temperature T even in a polydisperse case. I n a . fo rced inelastic system the repartition 
of energy seems to be very different and depends on the type of forcing usedtaEl 

In our system, the mean square velocity v 2 of particle i depends on its mass rrii and also on the proportion of all 
species j and their masses rrij . In all cases v 2 is constant at large time. Energy balance in a bidisperse system means 
that the agitation energy per unit time for a given species equals the energy lost in collisions with particles from all 
species. This can be written as follows: 



P(m s ,m s )w ss m s v 2 
P(mb,mb)wbbmbV b 2 



- P(m s ,mb)w sb m s v 2 
P(rrib, m s )wbsmbvl ■ 



mbVo, 



(4.1) 



where P(mi,mj) is the mean relative loss of kinetic energy by particle i when particles i and j collide. Clearly, this 
term depends both on mass and relative velocity. We also expect that P(m s ,m s ) should be equal to P(mb, mj). 
Therefore the quantities P(m s ,m s ) and P(mb,mb) calculated in a bidisperse system must be the same as those 
obtained for the monodisperse case. This term P(mj, m,) is therefore only a function of coefficients of restitution and 
friction. The term w^_represents the frequency of collisions of the grains i with j grains. We note that n s w s b = ribWbs- 
Using Enskog theoryEH, the frequencies of collisions of the grains in our two-dimensional system are given by: 




w S b = xV^(R 
Wbb = xV^iZRb 
_ Wbs = xV^{R. 

X is a cc*section factor and corresponds to the local radial distribution around a particleliii. 
shown inEJ that x does not depend on the type of particles but only on the packing fraction. 

The two limit cases, x s = and x s = 1, correspond to monodisperse situations with R 
respectively. In these two cases, we have determined numerically the four parameters v 2 (x s 



(4.2) 



We have previously 



R s and R 

0), v 2 s (x s 



Rb 

1). 



(x s — 0) and v 2 (x s — 1) by simulating a particle of radius Ri in a sea of particles of radius Rj. Using Eq. (fO|), 
we can calculate for these limiting values of x s the four parameters P(m,,mj). We have verified that P(nii,mi) is 
independent of the type of particle. We have found that P(m s ,m s ) = P(nib,mb) — 0.145, P(mb,m s ) — 0.229 (at 
x s = 1) and P(m s , m&) = 0.066 (at x s — 0). 

To a first approximation, we consider the P(mi,mj) to be independent of the relative velocity. We compare in 
figure uXl the values of v 2 and v 2 obtained from the numerical simulations and those deduced from Eq. (|4. 1[) . These 



values were calculated for different :r s at a packing fraction of 40%. The dashed lines correspond to the theoretical 
values and the symbols to the simulations. The agreement between simulations and theory is quite good, in particular 
for the big particles. 

However, the e nerg y lost in a collision does depend on the relative velocity, and therefore on v 2 and v 2 . To treat this 
correctly in Eq. (4.1), we show in figure [ll] that both, the kinetic energy of the system and v 2 /v 2 , decrease linearly 
with x s . Recalling that x s lies in the interval [0,1], we can thus write: 



Cx, 
;r./, ,J 



m s v 2 (x s ) + C{ l R * s) m b vl{x s ) = -^m h vl(x s = 0) + ^m s v 2 (x 



c 

it Hi 



m b v 2 (x s = 0) 



o) 



v 2 (x s ) _ V>(. 
vHxs) ~ v({x s =0) 



+ 



:(g.=o) 



(4.3) 



Using the four equations 4.3 and 4.1, we can calculate directly v 



'b ■ 



~i s ,mb) and P(rrib,m s ) as a function of 
x s - Note that different values of x s correspond to different values of the ratio v 2 /v 2 . The solid lines in figure [l(] 
correspond to the theoretical velocities squared calculated with this approach. The agreement with simulations is 
now very good. The values of P(m s ,mb) and P(mb,m s ) from the solution of our four equations are shown in figure 
|j~2| . The values of P(m Sl mb) for x s near 1 an d P(mf,,TO s ) for x s near should be taken with precaution, because 
their weight in the balance of energy (Eq. flj.lD ) is negligeable for x s ~ 1 and x s ~ 0. Indeed the energy loss in these 
limit cases corresponds to rare collisions. 

Note that the approximation of P(mb,m s ) by a constant (as done in the first approach) is fairly good. However, 
the value of P(m s ,mb) varies significantly with x s . This explains the slight difference between simulation and theory 
in the first approach. Finally we see that at large x s (corresponding to small ratio v 2 /v 2 ) the small particles gain 
energy in collisions with big particles. This phenomenon occurs only in dissipative forced gases where v 2 is no longer 
proportional to 1/m. 
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B. Non-uniform agitation 



We will now treat the case of non-uniform agitation. To do this we use exactly the same algorithm but introduce 
a gradient in the agitation by imposing the following spatial dependence for rj^: 

T)l(x) = a + bx, (4.4) 

where a and b are constants. The initial configuration of the system is taken to be the stationary state found in the case 
of homogeneous agitation. The gradient of the agitation is chosen such that the mean value of t/q (x) over the whole 
system corresponds to the agitation of the initial state. The simulations show that a concentration gradient appears in 
the system. The system reaches a non-uniform steady state where the density gradient remains present in the course 
of time. Note, however, that if we follow the motion of a particle, we find that it does visit the whole system. All 
the previous relations and equations are still valid in this case but one should consider a local "equilibrium" . Figure 
|l3| represents a typical configuration obtained at large time. The hot agitation is on the right side of the system. 
In the stationary state, the gradient of concentration balances the gradient of agitation such that the pressure is 
homogeneous throughout the system. 

The other main observation is that segregation appears in the system. The local proportion of s and b grains is 
no longer the initial one. The big particles are more sensitive to the gradient of agitation than the small ones. We 
show in figure |lj the local packing fraction for both species as a function of the position x in the system. We have 
found in all cases we have investigated so far that the stationary state always exhibits segregation. Our numerical 
results on segregation (the big particles are more concentrated ia the colder region) is in complete agreement with 
the theoretical calculations based on the granular kinetic theoryo. 



V. CONCLUSION 



The main purpose of this paper is to investigate granular mixtures numerically. We used an algorithm which keeps 
the dissipative particles agitated by applying to the grains an external random acceleration independent of their mass. 
In the case where the external acceleration is independent of the position of the particles (homogeneous agitation) we 
have shown that the system reaches a well mixed homogeneous stationary state even in the bidisperse case. We have 
therefore shown that agitation, and therefore diffusion, is an efficient mechanism for mixing. We have established 
a theoretical expression between the typical time of mixing, T m i%, and the number N Si b of collisions between s and 
b grains per unit of time which is valid for monodisperse as well as bidisperse systems. For the monodiperse case 
we have given the exact expression of T m i X as a function of the diffusion coefficient. For the bisperse case we have 
found that T m i x depends strongly on the initial configuration of the system. We have characterized the steady state 
reached by bidisperse assemblies and in particular we have established energetic relations which allow us to evaluate 
the square velocities of the grains as a function x s (the relative proportion of both species). 

In addition we have investigated the case where there is a gradient of agitation through the system and have shown 
that segregation appears. The main cause of this segregation is not related directly to the nature of the grains (size, 
dissipation, roughness) but originates from the presence of a temperature gradient. In a lot of mixing experiments 
(rotating drum, vibrated system) there often exists a gradient of agitation which then leads to a segregation process. 
Of course there are other types of segregation mechanisms which are sometimes more efficient, but this one is in some 
sense intrinsic. 
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FIG. 1. Configurations at three different times t = 1, t = 50 and t — 199. a): for an initially mixed system and b): for an 
initially segregated system. 

FIG. 2. Number of collisions per unit time between the two species as a function of time. The dashed (solid) line corresponds 
to the case of FIG. |l|a (b). 

FIG. 3. N 3t b(oo) — N Bt b(t) vs time. Solid line is the numerical average over 5 simulations. The dashed line is given by 
N St b{oo)exp{—t/T m ix), Tmix = 43.5. 

FIG. 4. Dependence of T mix on D. Q: value from simulations. Solid line: theoretical prediction. 

FIG. 5. Evolution, at short time, of an initially segregated bidisperse system; From left to right t = 1, t = 6 and t = 36. 
R a = 0.4 and; Rb = 0.6, the total packing fraction is 40%. x s = 0.5. 

FIG. 6. The local packing fraction of the small particles as a function of x (the lateral position) and t. 

FIG. 7. Vg(t) and Vb(t) as functions of time, for an initially segregated system. 



FIG. 8. N 3t b(t) vs t for a bidisperse gas. R a = 0.4, Rb = 0.6, and x a = 0.625. T m ix — 30 is obtained by fitting Eq. (3.8), 
shown as a dashed line. 



FIG. 9. r m ix for the bidisperse case as a function of x 3 . R 3 = 0A,Rb = 0.6 and C = 0.4. T m ix is constant. 

FIG. 10. The mean square velocities reached in t he stationary state as a function of x s . v 3 (□) and v 2 (O) obtained from 
simulations. Dashed lines are obtained from eq 4.1 assuming constant P(mi,mj)'s. Solid lines are the results from Eq. ( |l.l| ) 



and Eq. (4.3) 
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FIG. 11. a) The total kinetic energy of the system as a function of x s for R s = 0.4, R b = 0.6, C = 0.4 and ?jg = 22.0. b) 
vljvl as a function of x s for the same parameters. 



FIG. 12. Values of P(m s ,irib) and P(mb,m s ) vs x s calculated using eqs 4.3 and 4.1 for the same mixtures as in FIG. |ll| 
Dashed line: P(mb,m 3 ). Full line: P(m 3 ,mb). 

FIG. 13. Typical configuration obtained in the presence of a gradient in the granular temperature. R 3 = 0.4, Rb = 0.6, 
x a = 0.5, a = 4 and b = 0.9. C = 0.3. 

FIG. 14. Local packing fraction as a function of x for s (□) and b (O) grains. The mixture corresponds to FIG. [l3|. Averaged 
over the whole system, x s = 0.5 and C = 0.3. 
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Figure l.a 
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Figure l.b 
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